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FILTERING  BY  USE  OF  LOCAL  STATISTICS 


SOME  RESULTS  IN  NC1SE 

P.  F.  Yan* 

P.  F.  Li** 

Department  of  Electrical  Engineering  w 
Southeastern  Massachusetts  University 
North  Dartmouth,  Massachusetts  02?47 

In  a  recent  paper  jlj  J.  S.  Lee  proposed  a  method  of  using 
local  statistics  in  noise  filtering.  The  basic  assumption  is  that 
the  auffimble  mean  and  variance  of  a  pixel  is  equal  to  the  local 
mean  and  variance  of  all  pixels  within  a  fixed  range  surrounding 
it.  So  the  gray  level  of  any  point  of  the  original  image  can  be 
estimated  by  the  local  statistics  and  the  noisy  gray  level  of  that 
point.  This  implies  that  the  image  is  locally  stationary  and 
ergodic.  Though  this  assumption  is  not  always  exact,  it  has  been 
shown  to  be  quite  feasible  by  some  experimental  results.  We  have 
implemented  this  method  (for  simplicity  we'll  call  it  L-method 
hereafter)  in  our  minicomputer  (PDF  -  11/45). 

1.  Determination  of  Variance  of  Noise 

In  order  to  calculate  the  local  statistics  of  the  original 
image  the  variance  of  the  additive  noise  G^“  must  be  known. 

Sometimes  only  the  noisy  image  ia  available.  In  the  case  where 
the  image  is  composed  of  objects  and  a  nearly  uniform  background. 

0^  can  be  calculated  from  that  window  which  has  the  least  mean 

value.  As  the  mean  of  noire  is  zero,  the  window  with  the  least 

2 

mean  must  be  the  uniform  background  and  its  variance  must  be  . 

2 

In  the  experimental  study  an  uniform  noise  ( ^  =300)  was  added. 
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This  value  was 
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0^  calculated  by  above  mentioned  method  is  2?0. 
used  in  processing  the  noisy  image.  The  result  is  shown  in  Fig.  1. 
2.  Relation  to  Other  Filtering  Methods 

First  let  us  calculate  the  mean  square  error  of  the  L-method. 
Assume  that  the  signal  x  and  noise  n  are  independent  Gaussian 
random  variables,  n  has  zero  mean,  and  z  is  the  measured  gray  level. 
Then  we  have 


Z  •  .  =  x  •  :  +  n 
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The  conditional  probability  density  p(x/z)  must  be  Gaussian  with 
mean[x"]  and  variance  [O'*] : 
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where  x  is  the  mean  of  and  z=x  f 
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estimation  x  of  x- .  is  equal  to  [x j ,  this  is  the  formula  used  t 
estimate  x^  in  the  L-method.  The  mean  square  error  of  this 
estimator  E(\\X-i\\jz  )  is  simply  the  variance  of  p(x/z), 
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is  the  signal -to-noise  ratio. 


A.  Relation  to  Local  Averaging; 

In  the  local  averaging  method  we  replace  every  x^  by  x. 
the  average  square  error  ^  must  be 
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uniform  hence  x  is  the  optimal  estimator  of  x,  When  L-method 
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is  superior  to  local  averaging. 

B.  Relation  to  the  Kalman  Filter 

In  Kalman  filter  we  estimate  x^  by  recursive  estimation.  For 
the  points  that  do  not  have  neighborhood  points  information  the 
filtering  formula  becomes  ^See  ! j\  Chapter  ?  formula  (142,143). 
Also  [ 4 1  formula  (7*46))  x.  .=  x  +  k(z-x) 

'  S  1  J 


k  = 


Arc  < 


or* 


where  _  _  _  _ 

Kl{(o,o)  +  CrS  0**  r  CR2 

which  is  the  same  formula  used  in  the  L-method.  (See  equal, ion ( 1) ) . 

We  can  conclude  that  L-method  is  nothing  more  than  a  special  case 
of  the  Kalman  filtering  and  its  effect  is  inferior  to  Kalman  filter 
in  the  global  operation.  However,  its  computing  algorithm  is 
simpler  and  can  be  directly  implemented  for  real-time  processing 
where  a  parellel  processor  is  used. 

The  effect  of  filtering  is  dependent  upon  the  characteristics 
of  the  images  under  processing.  If  the  image  is  nearly  uniform  or 
the  gray  level  of  one  pixel  is  independent  of  other  pixels,  then 
local  mean  is  the  best  estimator.  If  the  gray  level  of  one  pixel 
is  related  to  other  pixels  but  the  relationship  is  unknown,  the  L- 
method  is  the  best  estimator.  If  this  relationship  can  be  assumed 
as  a  suitable  model,  then  Kalman  filter  is  the  best  estimator.  If 
the  image  has  more  complex  characteristics,  we  must  use  some  non¬ 
linear  methods.  Here  "best”  means  the  least  mean  square  error. 

3.  An_ Alternative  Approach  in  the  Case  of  Multiplicative  Nois_e 

Under  multiplicative  noise  the  degraded  pi/el  can  be  represented 

by  z . . =  x . .  u . .  ( 3 ) 

J  ij  ij  ij 

where  u. .  is  noise.  In  Ref.  1,  this  problem  is  treated  by  a  linear 
1  J 

relation  between  z,  x  and  u: 
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z  '  -  -  -  Ax,  .  +  Bu .  .  +  C 
J  ^ 

A,  B  and  C  are  determined  such  that  the  mean  square  error  between 


z  —  and  zl.  is  minimized.  The  resulting  Tillering  formula 


is 


X  •  •  =  X  +  K  •  •  ( Z  •  •  -  a  Y.  )  ,  .  ,  .  -  , 

lj  lj  lj  wmch  is  x  plus  some  correction  terra. 

Intuitively  for  the  multiplicative  noise  Inis  formula  is  not 

reasonable  as  in  the  additive  noise  case.  Therefore  we  propose 

.  .  A 

another  approach.  Intuitively  the  reasonable  form  ol'  x  may  he 


x .  .  =  lYZ  .  . 

iJ  ij 


Now  we  take  log  of  (3)t 


and  assume  lh  x  be  a  linear  function  of  In  z  i.e.  In  x  =  c  In  z  +  k 
(or  x  =  e  z  )  .  Here  c,  k  are  determined  such  that  the  mean  square 
error  between  In  x  and  lrfx  ,  E(ln  x  -  lrPx)2  is  minimized. 


C  -  £  -  (tin  X)2 

E  InZ  -  (Blnz)* 


t  In*  -  ( £  I^X) 
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The  problem  is  to  calculate  Eln  x,  Eln  x...etc.  li‘  the  statistics 

__  *2  _  1. 

of  noise  is  known,  i.e.  we  know  u  and  O^then  the  x  and  0^  can  be 
determined  from  the  following  equations: 


of.  _Sli 


l nee  Appendix) 


Here  z  and  CT*  are  the  local  mean  and  the  variance  of  the  noisy 
image . 

With  a  reasonable  assumption  of  the  probability  density  of  x, 

2 

Eln  x  and  Eln  x  can  be  calculated  from  £!nX  = J  in 

t  In  X  •=  jin  X  Jf>(x)dx 

*It  is  the  same  principle  as  Homomorphic  filtering.  Here  the 
processing  is  in  the  spatial  domain  using  the  local  statistics. 


Here  we  simply  list  the  resulting  filtering  formula.  And  1'or 
detail,  please  see  the  appendix. 

$■=  FinX  =  y  -  H  -/.S’ 

r=  tU\  =  k  -±f  *H(3-lnA)  +  3. S' 
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Experimental  results: 

The  multiplicative  noise  is  uniform  distritution  from  0.2  -  1. 
After  processing,  the  result  is  shown  in  Fig.  2A.  For  the  case  of 
combined  multiplicative  and  additive  noise,  the  result  of  processing 
is  shown  in  Fig.  3-  For  comparison,  with  the  same  noise, the  result 


where 
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of  using  the  L-method  is  shown  in  Fig.  2B. 
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1,  Derivation  of 

o <%*=Fxx-C£x)  ,  <K  -  t  u  -  (fu ) 
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2.  Calculation  of  Elnx,  Eln^x 

In  practice  the  distribution  is  usually  assumed  to  be  Gaussian. 

For  simplicity  we  consider  a  triangular  distribution  (Fig.  4)  as  an 

approximation  to  the  Gaussian  distribution.  The  distribution 

function  of  x: 
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By  using  the  formulae, 

J~ln  X  d  X  —  X  In  A  —  X  +  C  .  Jjf/nxdx  =**  L^-^rJ+C 
)xlnlx  dx  ~^lnlx  -  x'l  ~  J  ]  /  jln\4x  -  Xin\~  2x)nx  ->  2-r 

the  results  can  be  easily  derived. 

B.  Calculation  of  Illnu,  bln‘u 

When  the  noise  u  is  uniformly  distributed  from  a  to  b,  then 
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F  i  g .  1 

Filtering  of  noisy  image 
with  additive  noise,  G^*-»3oo, 
oi’  uniform  distribution 
(-30,  3C).  The  upper  figure 
in  the  intensity  profile 
along  a  scan  line  of  the 
image  in  the  lower  figure. 
Original  image  in  the  left, 
noisy  image  in  the  middle, 
filtered  image  in  the  right. 

Image  is  from  USC  data  base, 
//USC-025co-15. 


Fig.  2 A 

Filtering-  of  r.oiry  image 
with  multiplicative  noise 
of  uniform  distribution 
(0.2-i.)  by  using  our 
rnetliod . 
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Fig.  2B 

Filtering  of  noisy  image 
with  multiplicative  noise 
of  uniform  distribution 
(0.2-1.)  by  using  L-method. 


The  improvement  of  our 
method  over  the  L-method 
is  evident  by  comparing 
Fig.  2 A  with  Fig.  2B. 
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